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A conceptually appealing and computationally economical course-grained molecular-orbital (MO) theory for 
extended quasi-linear molecular heterostructures is presented. The formalism, which is based on a straightfor- 
ward adaptation, by including explicitly the vacuum, of the envelope-function approximation widely employed 
in solid-state physics, leads to a mapping of the three-dimensional single-particle eigenvalue equations into 
simple one-dimensional hole and electron Schrodinger-like equations with piecewise-constant effective poten- 
tials and masses. The eigenfunctions of these equations are envelope MO's in which the short-wavelength 
oscillations present in the full MO's, associated with the atomistic details of the molecular potential, are 
smoothed out automatically. The approach is illustrated by calculating the envelope MO's of high- lying oc- 
cupied and low-lying virtual ir states in prototypical nanometric heterostructures constituted by oligomers of 
polyacetylene and polydiacetylene. Comparison with atomistic electronic-structure calculations reveals that 
the envelope-MO energies agree very well with the energies of the ir MO's and that the envelope MO's describe 
precisely the long-wavelength variations of the it MO's. This envelope MO theory, which is generalizable to 
extended systems of any dimensionality, is seen to provide a useful tool for the qualitative interpretation 
and quantitative prediction of the single-particle quantum states in mesoscopic molecular structures and the 
design of nanometric molecular devices with tailored energy levels and wavefunctions. 
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I. INTRODUCTION 

Extended molecular structures appear in condensed- 
matter physics, macro- and supra-molecular chemistry, 
molecular biology, materials science, molecular electron- 
ics and optoelectronics, and nanotechnology. Rang- 
ing from one-dimensional to three-dimensional, these in- 
clude, among ma ny others, inorganic, organic and biolog- 
ical polymer and wires'^ dendrimersp' surfaces^ and 
graphene-based structures^ fullerenes^ and nanotubesP 
clusters^ and crystals \ n addition, there exist as- 
semblages of different fragments or materials, like hybrid 
structured' and heterostructures.^ 

When such an extended system is sufficiently large so 
that, for practical purposes, it can be considered as infi- 
nite, and periodic, so that a repeat unit cell can be iden- 
tified, periodic boundary conditions can be employed, ef- 
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fectively reducing the calculation of its electronic struc- 
ture to the one of the unit cell at a set of sampling k 
points. Within (Schrodinger or Kohn-Sham) orbital ap- 
proaches, this leads to a description in ter ms of energy 
bands and (Bloch) crystalline orbitals!- 1 * 5 * 10 * On the other 
hand, if the system lacks translational symmetry alter- 
native approaches must be sought. 

Non-periodic extended molecular structures of interest 
in the abovementioned fields can possess from hundreds 
up to millions of atoms. Depending on the degree of accu- 
racy required and the computational resources available, 
there are semi-empiricaf^ and ab-initio^ electronic- 
structure methods to choose from. When coupled 
with linear-scaling algorithms,^ these methods can han- 
dle large numbers of atomsp^ like the impressive re- 
cent million-atom implementation of orbital-free DFTP^ 
Other approaches, like density-matrix renormalization 
theory,^ elongation methods,^ fragment molecular- 
orbital theorypSland embedded cluster methods,^ which 
can deal, in principle, with large non-periodic systems 
have been formulated. In spite of all these developments, 
the calculation of the electronic structures of such sys- 
tems remains computationally demanding, which makes 
the formulation of accessible approaches for this task a 
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continuing effort in quantum chemistry. 

Naturally, what makes such electronic-structure meth- 
ods for non-periodic extended systems computationally 
expensive is their processing of information at the atomic 
(sub-nanometric) length scale. However, as the meso- 
scopic scale is approached, it can be expected that some 
physical (e.g. electronic and optical) properties of the 
system will depend on information at a larger (nano- 
metric) scale. This means that these properties could 
be obtained from a course-grained average of the atom- 
istic information. Hence, it is of considerable interest to 
devise methodologies that permit the direct determina- 
tion, with good accuracy of such properties bypassing 
the atomistic details. This would amount not only to 
large computational savings but also to increased under- 
standing, since confounding atomistic information would 
be automatically averaged out. 

An approach of this guise was introduced early in 
quantum chemistry for the determination of the elec- 
tronic structure and spectra of lineal!^ and cyclical con- 
jugated molecules, namely the 'free-electron molecular- 
orbital' (FEMO) model. This approach, and exten- 
sions thereof,^ were applied to small- and medium-size 
molecules. FEMO models have proven helpful as well 
in the study of other molecular properties, for instance 
the correlation of photoionization resonances with bond 
lengths^ and the densities of states in fullerenesP^ 

An approach of the same spirit was also independently 
developed early in solid-state physics for the treatment of 
impurities and other slowly-varying weak perturbations 
in bulk semiconductors, namely the 'envelope-function 
approximation' ( EFA) , also called 'effective-mass approx- 
imation' (EMA) .12321 Later, this method was success- 
fully, albeit heuristically, applied to semiconductor quan- 
tum wells embedded in a bulk material.^ Since then, 
the EFA has provided a very useful framework for the 
interpretation and prediction of the electronic, optical 
and conduction prop erties of heterostructures built from 
crystalline materials l 10 l 12 l 3 ° l due to its relative concep- 
tual simplicity and economical computational implemen- 
tation in com parison with ato mistic electronic-structure 
methodologiesPS ^ 14 ' 16 ' 17 ' 30 ' 

Although the EFA and the FEMO model have a lot in 
common, in particular the use of effective box poten tials, 
the first has been put on a hrm theoretical basis 
whereas the second is based on (very clever) ad hoc con- 
siderations introduced on largely intuitive grounds. 

The main goal of this paper is to demonstrate that 
the formalism and implementation of the EFA can be 
straightforwardly, but rigorously adapted for the deter- 
mination of the course-grained electronic structure of ex- 
tended, albeit hnite, molecular systems. This approach is 
hereafter called 'envelope molecular-orbital (EMO) the- 
ory'. Attention will be focused on organic quasi- linear 
heterostructures, for which it will become apparent that 
the EMO method can readily handle systems with meso- 
scopic lengths. 

Section [H] begins with a compact review of the EFA 



formalism in the context of a polymer heterojunction, 
put in a language more quantum-chemical than found 
in the solid-state physics literature. The authors do 
not claim originality as to the formulation, which is 
standard, but hope that their exposition will be help- 
ful both to chemical physicists/quantum chemists, who 
may be unfamiliar with the subject, and solid-state physi- 
cists interested in molecular nanostructures. Next, it is 
shown how the EFA can be rigorously adapted for (finite) 
molecules (EMO theory), by including explicitly the vac- 
uum in the formalism. In Section |III| illustrative proto- 
types of quantum dots (QD's) are designed by combining 
oligomers of irans-polyacetylene (PA) and polydiacety- 
lene (PDA). The EMO input parameters required for 
such design are extracted from electronic-bandstructure 



calculations of these polymers in Section IV The results 



of the EMO calculations for the prototypical QD's are 
presented, discussed and compared with the results of 
atomistic electronic-structure calculations in Section IV] 



Finally, Section VI provides concluding remarks and per- 
spectives for future developments and applications of the 
EMO theory. 



II. ENVELOPE MOLECULAR-ORBITAL THEORY 

As a prelude, let us consider a 7r-conjugated polymer 
A constituted of an infinite quasi-linear chain of iden- 
tical monomers a of length l a aligned in the z direction 
(Fig. [TJa)). As in any periodic system, the single-particle 
states, or crystalline orbitals (CO's), of such polymer 
are organized into continuous bands.^' This work fo- 
cuses on polymers with the highest occupied (jr) band, 
or valence band (VB), and the lowest unoccupied (it*) 
band, or conduction band (CB), separated by an energy 
gap, i.e. semiconductors or insulators. For example, the 
bandstructures of PA and PDA are displayed in Fig. [2] 
According to Bloch's theorem, each CO can be written 
asSE 



ip n {k,r) = e lkz u n (k,r), 



(1) 



where n is the band index, k = k z is the wavenumber 
conjugate to the longitudinal direction z, and u n (k,r) 
is the background Bloch function, which possesses the 
same spatial periodicity as the molecular potential. Since 
ijj n (k,r) is a periodic function of k with period 2tt / '£ a , 
this variable is considered to be restricted to the first 
Brillouin zone (BZ), |fc| < n/e a . The CO's of polymers 
are commonly calculated by quantum-chemical methods, 
employing linear combinations of basis Bloch orbitals, 
which, in turn, are translational-symmetry-a dapted lin- 
ear combinations of atomic orbitals (LCAO's)! 1 * 5 * 14 ! 

Now, let us consider a 7r-conjugated heteropolymer A- 
B formed by joining two semi-infinite strands of parent 
homopolymers A and B (Fig. [ljb)). Due to the break- 
ing of the spatial periodicity at the junction, Bloch's 
theorem (Eq. ([!])) no longer applies and, consequently, 
the band structure of single-particle states is destroyed. 
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FIG. 1. (a) The infinite homopolymer A. (b) The infinite 
heteropolymer A-B. (c) The finite block co-oligomer A-B-C, 
including the lateral vacua. 
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FIG. 2. EH bandstructures of the parent polymers. Due to 
k/ — k symmetry only half of the Brillouin zone is shown. 



Thus, these states must be described in terms of molec- 
ular orbitals (MO's) instead of CO's. In a conventional 
quantum-chemical approach, the MO's are represented as 
LCAO's, or, if this system happens to belong to one of the 
point-symmetry groups, as linear combinations of basis 
symmetry orbita ls, wh ich, in turn, are point-symmetry- 
adapted LCAO'sPHS (In practice, since this hypotheti- 
cal heteropolymer is infinite, artificial asymptotic or pe- 
riodic boundary conditions must be imposed.) 

Alternative approaches can be formulated, in which 
the MO's of the A-B heteropolymer are constructed in 
terms of the CO's of the A and B parent homopolymers. 
For example, the Green's matrix formalism (GMF)J22lhas 
been developed for the atomistic determination of local- 
ized interface states. The EMO theory is in the same 
spirit, but it yields a non-atomistic (course-grained) ap- 
proximation to the MO's. 

The EMO formalism begins with a spatial partitioning 
of the single-particle Hamiltonian, as follows. Let us de- 
note by Ra and Rb the spatial regions occupied by the 
semi-infinite strands of A and B, and by zab the z co- 
ordinate of their junction, which can be regarded as the 



midpoint of the A-B link (Fig. [ijb)). In the R A region 
(z < zab) the Hamiltonian is expressed as 



H&a) =H a Ra) +U {RA \r) 1 



(2) 



where H a Ra ^ is the single-particle Hamiltonian of the 
parent homopolymer A and U^ Ra> {t) is a perturbing po- 
tential defined so as to take into account the changes in 
the interactions acting on the particle in region Ra due 
to the replacement of the semi-infinite strand of A by the 
semi- infinite strand of B in the Rb region (z > zab)- By 
definition, this potential is identically zero for z > zab- 
For example, if Ha is the closed-shell Fock operator,^ 
H A = f A = T + 

VA.core 

(r) + v A ,cx(r), where the second 
and third terms are the core and Coulomb-plus-exchange 
potentials, then 



U( RA \r) = [v 



(Rb) 



vfex^] [1-Q{z-zab)], 



(Re 



(r) 



(3) 



where 0(z — zab) is the Heaviside unit step function.^ 
Although, for the present purposes, the precise form of 
the perturbing potential need not be specified, it is evi- 
dent that its magnitude is a maximum at the junction 
and usually small everywhere in comparison with the 
magnitude of the effective potential present in H a Ra ^ (in 



the example, v 



(Ra 
A, core 



Pj(r)). Moreover, since the 



7r electron-nucleus and electron-electron interactions are 
screened, this potential can be expected to be of relatively 
short range. When a partitioning analogous to Eq. ^ is 
defined for the Rb region, the full single-particle Hamil- 
tonian can be written in the spatially-piecewise fashion 

H = [1 - Q(z - zab)} + H^ Rb) Q{z - zab)- (4) 

The CO's of the parent polymers J = A,B are the 
eigenfunctions of respective single-particle Hamiltonians, 



Hj^{k,f)=e^{k)^ n J \k,f), 



(5) 



and, consequently, constitute a complete orthonormal 
set, 



R* 



d\^>{k',r)^ n J \k,r) = 5{k' - k)S mn . (6) 



Therefore, the MO's of the A-B heteropolymer, which 
are the eigenfunctions of the single-particle Hamiltonian 

0. 



(7) 



can be represented in each region in terms of the CO's of 
the corresponding parent homopolymer. In the Rj region 
such 'linear combination of crystalline orbitals' (LCCO) 
reads 



(8) 



(f) = E / dk fn Rj) ( fc > 

= W dke^fi R -'\k)uii J \k,f), 

„ J BZ 
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where the amplitudes f n Rj \k) are assumed to possess 
the same fc-periodicity as the CO's and Bloch's theorem 
(Eq. 0) was used in the second equation. Substitution 
of this expansion into Eq. ^ followed by projection onto 

iftm (k',f) and Fourier transformation to z-space yield 
the basic EMO equations 



e^{K)F^){z) + Y J I dz'U^{z,z')Fi R -'\z') 

n J R 



(9) 



erator K — —i4- in the dispersion relation of Eq. ab\, 



U {Rj) (z z') = -i- 

u ran V'i ) — ^ 



dk 



dk'e lk ' z e- lkz ' 



BZ J BZ 



(10) 



and F^ n Rj \z) is an 'envelope MO' defined by 



27r -IBZ 



dke lkz f£>Hk). (11) 



with j = a, b. (It should be noticed that the functions 
\$>Rj (jfj anc j pRj ( z ) are defined over the entire z— space, 
the super-index indicating that only their values confined 
within the Rj region are being considered.) In obtaining 
the first term of Eq. Q it was recognized that 



where Eqs. Q, ([2]), ^ and ^ were taken into ac- 
count, im\k) is a 'generalized kinetic-energy operator' 
obtained by formally replacing the variable k by the op- 



2tt 



BZ 



dke ik *eg(k)f™(k) = / dz'i^(z')F^(z-z') 



dz'e^{z')t 



-iz'K 



i^(k)Fi R -'Hz) 7 



(12) 



where in the first line the convolution theoreirPfil was 
used, with 



2tt 



dke tkz e^{k), 



BZ 



the second line the translation operator^ 

(z - z>) 



(Rj) 



f(z')F^ J) (z) = cM-^'K)F^ J \z) 
was introduced, and in the third line the inverse of Eq. 
(131 was used formally replacing the variable k by the 
operator K = —i4^. 

The set of real-space integro-diffcrential equations ([9| 
for the envelope MO's is equivalent to the Schrodinger 
equation Q, but much more difficult to solve in practice. 
Sophisticated techniques for dealing with the nonlocal 
terms [/^''(z, z') have been developed.^ Fortunately, 
these equations are also amenable to approximations that 
simplify drastically the determination of the MO's, at the 
same time retaining the essential features necessary for 
describing the physics of the system. Since the main 
purpose of this paper is to show that Eqs. ([9| can be 
fruitfully adapted to molecules, the second, much simpler 
approach is adopted here. 



The first and key approximation, which will be 
called 'uncoupled approximation', consists of neglect- 
ing altogether the nonlocal perturbation matrix ele- 
(13) ments ulnn\z, z'). This approximation is valid suffi- 
ciently far from the heterojunction, where U tyRj \r) is 
not large enough to couple the bands appreciably. The 
resulting uncoupled equations for the envelope MO's are 
free-particle-like Schrodinger equations with generalized 
kinetic-energy operators Em\k). Now, in both regions 
each MO (Eq. [8]) contains contributions from one band 
only, for which reason this approximation is also termed 
'one-band approximation'. It is reasonable to expect that 
the interfacial region around the junction, where this ap- 
proximation fails, has a width wab ~ (-a + &b- 

For the sake of concreteness, from now on interest will 
be focused on occupied (hole) MO's with energies near 
the VB edges (or highest occupied crystalline orbitals 
(HOCO's)), ei J) (k = k v ), and virtual (electron) MO's 
with energies near the CB edges (or lowest unoccupied 

crystalline orbitals (LUCO's)), ei (k — k c ), of the par- 
ent polymers, since their optical and conduction prop- 
erties are largely dominated by these ir and 7r* states, 
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respectively.^ Furthermore, let us assume for the time 
being that the HOCO and LUCO of homopolymer A 
lie higher and lower, respectively, than the HOCO and 
LUCO of homopolymer B (see Fig. Q). 

In general, three types of MO are possible in such A- 

( B) 

B heteropolymer: (1) For holes with e„ (k v ) < E < 
£y A \k v ) or electrons with ei A \k c ) < E < e{ B \k c ), 
MO's confined mainly to the (classically-allowed) Ra re- 
gion and decaying into the (classically-forbidden) Rb re- 

( B) 

gion; (2) for holes with E < e v (k v ) or electrons with 

( B) 

E > e\ >(k c ), MO's delocalizcd over the entire structure; 
and (3) MO's localized around the junction ("interface" 
states), which, evidently, cannot be treated within the 
uncoupled approximation and will not be considered in 
this work. MO's confined to the Rb region and decaying 
into the Ra region are impossible with the band align- 
ment presently assumed. 

Let us consider a type (1) MO with energy E. In the 
Ra region, and outside the interface, such MO must be 
very similar to the CO of parent homopolymer A with 
energy E — el (fc = fc B ), which implies that the func- 
tion fn B ' A \k) must be highly peaked around kE- Hence, 
since uh^{k, r) is a slowly-varying function of k, it is a 
good approximation to neglect its k dependence within 
the range of fn (k) and take it out of the integral in 
Eq. ([8]). On the other hand, in the Rb region at energy 
E there are not any states of the parent homopolymer B, 
which implies that the MO must comprise contributions 
from CO's of B with energies spanning, in principle, the 
entire BZ, in order to achieve the destructive interfer- 
ences necessary for the formation of an evanescent wave. 
Thus, the function fn RB \k) must be much wider around 
ks than fr R ' A \k) and the neglect of the k dependence of 
u\ B \k,r) in Eq. (JsJ) is unjustified. Nevertheless, since 
the evanescent portion of the MO is of short range, the 
relative error introduced by this approximation in the 
Rb region should be actually small. By the same to- 
ken, outside the interface a MO of type (2) with energy 
E resembles the CO's of parent homopolymers A and B 
with energies E — el^^fc^) and E = e^iks) in the 
Ra and Rb regions, respectively. In conclusion, taking 
into account Eq. |TT| ) and the assumed proximity to the 
respective band edges, in both regions hole and electron 
MO's of types (1) and (2) are represented as 



S(k — ks), where Fn RA \z) becomes a plane wave and 
Bloch's theorem ([lj is recovered in Eq. pi): since 
\k ma x\ = 7r/4 = W-Vnim then A mi „ = 24, indicating 
that the EMO's vary appreciably only over a range of sev- 
eral monomers, in contrast with ui^\k,r) which varies 
widely from one monomer to the next. Thus, Eqs. (14) 



and ( 15 ) indicate that an MO is approximately consti- 



tuted by a (rapidly-varying) background Bloch function 
modulated by a slowly- varying EMO, which is analogous 
to Bloch's theorem for a CO (Eq. ([I])) where a back- 
ground Bloch function is modulated by a slowly-varying 
plane wave. 

Let us now examine more closely the nature of the 
operator e { n J \k). Near the HOCO and LUCO, k„ 
(n = v,c), where the energy dispersion is usually nearly 
parabolic (see Fig. [2]), this function can be expressed as 

In or- 



in \k k n ) — £ n ^ (k — fc n ) + ^ ^ 



(k-k n y 



k=k n 



der to make an analogy with the free-particle dispersion, 
the 'effective mass' 



1 



d z e 



dk 2 



(16) 



k=k„ 



is introduced. (It should be noticed that, since the VB 
dispersion has the shape of an inverted parabola around 
the HOCO, the effective mass of the holes is negative.) 
Formal replacement of the variable k—k n by the operator 
K = —i£ yields 



i ( n J \k)^e n J \k = k n ) 



h 2 K 2 



2m 
h 2 



(j)* 



Ira 



(j)» dz 2 



(17) 



In this expression, k n — or k n — Tr/£ a , depen ding on 
how the VB and CB of the J homopolymer runPEl F or 
example, in Fig. [2]it can be seen that the frontier CO's of 
PA are located at the edge of the BZ, whereas the frontier 
CO's of PDA are located at the center of the BZ. Eq. (fl7|) 



constitutes the third approximation of the development, 

the familiar 'parabolic EMA'j 10 ' 12 ' 3 D whence i n J \k) be- 
comes the standard kinetic-energy operator, with the en- 
ergy measured from the band extremum. Now, Eq. ^ 
reduces to the free-particle Schrodinger equation 



^ h Rj \^=Fl R '\z)ui J \k Vl r), (14) 
¥ e Rj \r) S F^ Rj \z)u[ J \k c ,r). (15) 

These expressions constitute the second approximation 
of the development and will be termed "single-envelope 
MO approximation" . 

At this point, it is important to emphasize that, since 
fn R ' A \k) is highly peaked around k n , according to Eq. 
|llj) Fi RA \z) must be a slowly- varying function of z. 
This can be seen by considering the limit f n R (k) —> 



i J \k n ))F 7 { R 'Hz)^EF^\z). 



2ml J) * dz2 

(18) 

When both regions are considered together the follow- 
ing picture emerges: To the left and right of the inter- 
face the carriers behave like free particles with effective 
and m^* subjected to constant potentials 



masses m 



(A)* 



Sn (k n ) and e n (k n ), respectively. However, within the 



interface, where Eqs. (14), (15) and (18) do not hold, 



the behavior of these quasi-particles remains unknown. 
In order to obtain a well-posed Sturm-Liouville problem 
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for the EMO's across the entire structure, Eq. (18) is 
artificially extended into the interface. This fourth ap- 
proximation is reasonable because, as discussed above, 
EMO's of types (1) and (2) are expected to vary lit- 
tle across the interfacial region. This extension can be 
accomplished by introducing any smooth interpolations 
across wab of the band-extremum potentials and effec- 
tive masses between the pairs of values (k n ) ,s^(k n ) 



(A)* (B) 

and m n ,m„ , respectively. Now, Eqs. (18 1 for both 



regions can be comprised into the single equation valid 
across the entire structure 



1 



2 dz \^m* (z) dz 



+ V n (z) 



F n {z)^EF n {z), (19) 



where m* n {z) and V n {z) are the position-dependent ef- 
fective mass and band-extremum potential, respectively, 
and the BenDaniel-Duke expression for the kinetic- 
energy operator with a position-dependent mass^l has 
been used. 

This procedure automatically enforces the conti- 
nuity of the EMO's and their derivatives at the 
heterojunction!^ However, it should be mentioned that, 
according to Eqs. (14) and (15), the continuity of an 



one to each vacuum), since the others are too far away. 
The full single-particle Hamiltonian now assumes the 
spatially-piecewise form 

H = H(R«)[l-e{z-z aA )} (20) 

Y 

+ H (Ri) [&(z - - Q(z - z M+1 )] 

i=A 

Thus, the uncoupled approximation can still be used 
for all the oligomer-oligomer heterojunctions. On the 
other hand, this approximation breaks down at the 
vacuum-oligomer interfaces, for the following reason: 
The single-particle Hamiltonian in a vacuum region of 
the heterostructure, say R a , has the form H^ Ra > = 
f + {¥). In this case the magnitude of the per- 

turbing potential cannot be considered either small or 
short-ranged, since the effective potential in the Hamil- 
tonian of the "polymer" a is null. For example, if 
the effective potential of the A oligomer is the closed- 
shell Coulomb-plus-exchange potential, then U^ Ra '{f) — 



(r) + v A R cx(r) [1 — Q(z — Zqa)]- By the same 
token, the contribution of a to the perturbing potential in 



(Ra 
Alcove V 



EMO apparently does not guarantee the continuity of the Ra region would be 



the full MO at the heterojunction, since U„ (fc„,r) and 

Un^ (k n , t) do not, in general, match at that point .^Nev- 
ertheless, it can be shown that this is not the case, by 
taking into account the phase invariance of the Bloch 
orbitals (Eq. [IJP 

Finally, let us consider a finite n-block 7r-conjugated 
co-oligomer A-B- ■ ■ -Z formed by joining a series of n 
oligomers of parent polymers A, B, ■ ■ ■ Z. For exam- 
ple, Fig. [ljc) illustrates a three-block structure A-B-C. 
In general, three types of bound MO's are possible in 
this system: (1) MO's confined mainly to one or several 
oligomers and decaying into the other regions; (2) MO's 
delocalized over the entire structure; (3) MO's localized 
around a junction. Specific examples of type (1) and (2) 
MO's are found in Section IVl 

The EMO formalism for a single heterojunction can be 
extended to such molecules by making the following con- 
siderations: (1) The key idea is to regard the co-oligomer 
as an a/A/B/ ■ ■ ■ /Y/Z//3 heterostructure, where a and 
j3 stand for the semi-infinite left and right vacua, respec- 
tively (Fig. [l|c)). (2) To accommodate these vacua in the 
development, they are considered as hypothetical parent 
polymers endowed with empty lattices. (3) Now, the left 
and right chain ends can be formally regarded as a — A 
and Z — f3 heterojunctions, respectively. (4) Finally, each 
heterojunction is treated as if it were isolated from the 
rest. 

Assumption (4) is justified on the following grounds. 
The perturbing potential in a given region (oligomer or 
vacuum) now contains contributions from the effective 
single-particle potentials of all the other regions. How- 
ever, only the nearest-neighboring regions contribute ap- 
preciably to this potential (two to each oligomer and 



A, core 

Zoa), which is not small either. Nevertheless, since the 
a-A and Z-/3 interfaces are short in comparison with 
the length of the molecule, and an MO whose amplitude 
reaches a chain end must decay rapidly as it penetrates 
into the vacuum, this breakdown is expected to cause 
only small errors in practice. 

Consideration (2) means that each vacuum "polymer" 
possesses a parabolic (free-particle) dispersion relation 
around any k = k z and CO's with the forms tpn^ik, r) — 
e ikz u„ (x,y), where (x, y) = e lk * x e lk yV with 7 = oc,/3 
(see Eq. (fi])). Therefore, the single-envelope MO approx- 
imation (Eqs. (14) and (15)) and the parabolic EMA 
(Eq. (17)) are actually exact in the R a and Rp regions. 
Moreover, within the molecule the single-envelope ap- 
proximation continues to be valid outside the interfacial 
regions, and, for a given MO, the parabolic approxima- 
tion works best in regions where its energy is close to the 
corresponding parent polymer band extremum. 

Now, the hole and electron EMO's are determined 
across the entire molecule by Eq. (19), where the 
potential-energy functions and position-dependent effec- 
tive masses assume one-dimensional piecewise-constant 
profiles, except within the interfaces where smooth inter- 
polations of these functions are employed. 

Fig. [3] illustrates an example of a three-block co- 
oligomer of the A-B- A kind, together with the piecewise- 
constant intramolecular potential profiles for electrons 
and holes. An intramolecular barrier height for electrons 
or holes is given by the difference between the energies 
('offset') of the LUCO's or the HOCO's, respectively, of 
the neighboring parent polymers, 



T (PA) 



( k n 



e 



(PDA) 



{k n )\. 



(21) 
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drawn as 
shown. 



Structure and intramolecular electron and hole 
-energy profiles, according to the EH alignments, of 
gomer PDAg/PAr/PDAg. The hetero junctions are 
abrupt and the vacuum and LOCO levels are not 



ments of the energies of the LUCO's and the HOCO's, 
respectively, of PA and PDA. In Section [rv] they are ob- 
tained by means of two methods and the results are pre- 
sented in Tables [T| and [TT] and Fig. |4j On the other hand, 
the LOCO energies, required for handling the oligomer- 
vacuum junction for holes, are not available. However, in 
practice this does not pose a problem, since the energy of 
this state is so far down that its magnitude can be taken 
as infinite without incurring in any significant errors. 

Fig. [3] shows the resulting intramolecular profiles for 
a symmetrical three-block structure. For interpretation 
purposes, it is useful to notice that although the in- 
tramolecular profile for holes looks like a potential bar- 
rier, it can be thought of as a well if the effective mass of 
the holes is taken as positive. 

T hese qu asi-linear architectures constitute 
Qrj' jio | i2 | 30 | 45 | smce the- particles can be confined 
in the three directions. The design of molecules of 
this type can be useful for optical and electronic 
applications,^ since they can exhibit relatively localized 
states with controllable sizes and energies, as will be 
appreciated in Section [V] 

Studies of the electronic states of heterostructures 
similar to these have been reported by other authors, 
employing conventional quantum-chemical or solid-state 
methods!™! 



By the same token, a molecule- vacuum barrier height for 
electrons is given by the offset between the LUCO of 
the lateral parent polymer (in the example, PDA) and 
the LUCO of the vacuum (which lies at zero energy). 
On the other hand, a molecule-vacuum barrier height 
for holes cannot be defined in an analogous way because 
the vacuum does not possess a HOCO. Nevertheless, the 
molecule-vacuum junction for holes can be handled by 
taking into account that the hole spectrum of the lateral 
parent polymer is bounded from above and from below by 
the energies of its HOCO and lowest occupied crystalline 
orbital (LOCO), respectively, and that the hole effective 
mass is negative. The molecule- vacuum junctions are not 
indicated in Fig. [3j since the vacuum and LOCO levels 
lie very high up and deep down, respectively. Naturally, 
the magnitudes of the electron and hole effective masses 
in the vacuum regions are equal to the normal electron 
mass, mo. 



III. 



MODEL HETEROSTRUCTURES 



To illustrate the theory, simple prototypical n-block 
(n=l,2,3) structures were studied, employing the com- 
mon parent polymers PA and PDA. These are shown 
in the first columns of Tables |T] and |Hj where the sub- 
scripts indicate the numbers of monomers in the con- 
stituent oligomers and the vacua have not been indicated. 
The EMO intramolecular potential profiles for electrons 
and holes in these systems are determined by the align- 



IV. COMPUTATIONAL IMPLEMENTATION 

A practical implementation of the EMO theory re- 
quires four input parameters for each parent polymer: 
the HOCO and LUCO energies and the hole and elec- 
tron effective masses. The accuracy of the MO energies 
and envelopes produced by this approach depends on the 
accuracy of these parameters and the method of solution 
ofEq. (pi). 



Ideally, the input parameters are extracted from ex- 
perimental measurements. Specifically, according to 
Koopmans' theorem,^ approximations to the LUCO and 
HOCO can be obtained from photoelectron spectra. On 
the other hand, th e effectiv e masses can be obtained from 
the conductivities.^ 12 * 46 * In lieu of experimental data 
these parameters can also be determined computation- 
ally from the bandstructures of the parent polymers and 



Eq. (16), which is the strategy adopted here. 



The goal of this section and the next one is to demon- 
strate the internal consistency of the EMO approach. 
Hence, its predictions will be contrasted with the re- 
sults of atomistic electronic-structure calculations, not 
with experimental data. Moreover, for this comparison 
to be meaningful, such calculations are performed em- 
ploying the same method as the one used for determin- 
ing the bandstructures. Since the co-oligomers studied 
are rather long, in order to avoid unnecessary complica- 
tions the following simple methodology was employed. 

First the geometries of the PA and PDA monomers 
were pre-optimized employing the semiempirical Austin 
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TABLE I. EH values of the effective masses (in units of the vacuum electron mass, mo) and the frontier CO energies of the 
parent polymers, and the high-lying hole and low-lying electron energies of the n-block structures. A state marked with an 
asterisk is unbound with respect to the corresponding intramolecular well. SE = \Eeh — Eemo\- 
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(eV) 


(eV) 


(meV /mEh) 


PA 


0.041 


0.037 


HOCO 


-10.944 












LUCO 


-9.572 






PDA 


0.048 


0.034 


HOCO 


-11.338 
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-10.969 


1/0.037 
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-10.995 
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LUMO+2 


-9.359* 


-9.338* 


21/0.77 








HOMO-1 


-11.379* 


-11.412* 


33/1.21 


PDA20/PA20 






HOMO 


-11.143 


-11.152 


9/0.331 








LUMO 


-9.411 


-9.410 


1/0.0367 








LUMO+1 


-9.320* 


-9.248* 


72/2.65 








HOMO-2 


-11.203 


-11.475 


272/9.99 








HOMO-1 


-11.093 


-11.089 
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PA 65 






HOMO 


-11.020 


-10.980 


40/1.47 








LUMO 


-9.491 


-9.531 


40/1.47 








LUMO+1 


-9.416 


-9.409 
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LUMO+2 


-9.301 


-9.207 


91/3.34 



TABLE II. Analogous to Table I, employing DFT instead of EH. 



System m* h 


m* 


State 


DFT 


EMO 


SE 


(mo) 


(m ) 




(eV) 


(eV) 


(meV/mEh) 


PA 0.066 


0.063 


HOCO 
LUCO 


-4.120 
-3.207 






PDA 0.075 


0.065 


HOCO 
LUCO 


-4.771 
-3.740 










HOMO-2 


-4.656 


-4.834* 


180/6.6 






HOMO-1 


-4.515 


-4.708 


190/7.0 


PDA20/PA20 




HOMO 


-4.260 


-4.283 


23/0.85 






LUMO 


-3.582 


-3.689 


110/4.0 






LUMO+1 


-3.459 


-3.541 


82/3.0 






LUMO+2 


-3.317 


-3.305 


12/0.44 






LUMO+3 


-3.186* 


-3.057* 


129/4.7 






HOMO-2 


-4.982 


-5.247 
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HOMO-1 


-4.835 


-4.983 


148/5.44 


PDA20 




HOMO 


-4.735 


-4.824 


89/3.3 






LUMO 


-3.656 


-3.679 


23/0.84 






LUMO+1 


-3.541 


-3.496 


45/1.65 






LUMO+2 


-3.378 


-3.191 


187/6.87 
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FIG. 4. DFT and EH band alignments of the parent polymers. Dotted lines: LUCO's, solid lines: HOCO's. 



Model 1 (AMI), 47 as implemented in the Gaussian03 
package.BS This was accomplished by optimizing the ge- 
ometry of oligomers of increasing lengths, until all the 
bond lengths and angles of the central monomer con- 
verged. For simplicity, all the oligomers were forced to 
maintain a planar configuration. In all cases, convergence 
was well achieved with decamers!^ 

Second, taking these optimized monomers as unit cells, 
single-point bandstructure cal culatio ns were performed 
at the Extended Hiickel (EH f 1 50 1 51 1 and DFT^ levels. 
For the EH calculations the BICON-CEDiT packageP 
was employed, adjusting the value of the K parameter 
that appears in the bond integrals^ so that the calcu- 
lated bandgaps agreed with the reported experimental 
valuesP2 The optimal K values turned out to be similar 
for both polymers. Since these parent polymers are to be 
used for building heterostructures, a single value of K = 
2.43 was determined to provide a good compromise!^ 
The resulting bandstructures are displayed in Fig(2] For 
the DFT calculations the SIESTA packagtP 1 was em- 
ployed, using the PBE functional and the pseudopoten- 
tials and DZP basis set implemented in it. The resulting 
bandstructures are not shown. 

Third, the parabolic effective masses were extracted 
from these bandstructures, according to Eq. (16). In 



the EH case, it should be borne in mind that the ad- 
justment of K is expected to yield good bandgap values, 
but is not guaranteed to produce accurate bandwidths 
and, consequently, effective masses. Therefore, the lat- 
ter were determined more accurately by means of a real- 
space oligomeric extrapolation^ based solely on HOM O 
and LUMO data for oligomers of increasing sizesp^^ de- 
tails and examples of which will be provided elsewhere. 

Tables |I| and [H] contain the frontier CO energies and the 
effective masses of PA and PDA extracted from the EH 



and DFT bandstructures. In Fig. [4] it can be seen that 
the band alignments predicted by DFT are different from 
the ones predicted by EH. (The EH and DFT alignments 
obtained are co mmonly referred to as 'type F and 'type 
IF, respectivelyF 10 ' 12 ' 30 ! ) The absolute HOCO and LUCO 
values produced by DFT should be more accurate than 
the ones produced by EH, whereas the opposite should 
be true for the bandgaps, since EH was parameterized 
here to reproduce the experimental bandgaps and DFT 
is known to underestimate these values.^ Thus, at this 
point it cannot be reliably established which the correct 
alignments are. In addition, the effective masses pre- 
dicted by DFT are somewhat larger than the ones pre- 
dicted by EH, by at most a factor of two. In any case, 
what is important here is that now the performances of 
the EMO method employing parameters obtained with 
two very different levels of theory can be contrasted. 

Armed with the eight EMO input parameters, the po- 
tential profiles for electrons and holes of the prototyp- 
ical heterostructures were easily defined. In this con- 
nection, it is important to realize that the lengths of 
the constant segments in these profiles are calculated by 
adding the lengths in the z direction of the optimized 
monomer structures, not by adding the lengths of the 
bonds. Now, the hole and electron EMOs can be ob- 



tained by solving Eqs. ( 19 1 employing any of the meth- 
ods available} 5 -^ whose requirements of computer memory 
and time are practically negligible in comparison with 
standard atomistic electronic-structure methods. In this 
work, the purely numerical method described in refP^ 
was chosen. 

For the atomistic molecular electronic structures, 
single-point EH and DFT calculations were performed 
employing the implementations included in the Hyper- 
chem v.7.01^1 and SIESTA^H packa ges using the same 
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AMI pre-optimized monomer geometries, parameteriza- 
tion, functional and basis set as for the polymer band- 
structure calculations. 



V. RESULTS AND DISCUSSION 

A. EH Input Parameters 

In this subsection the EMO method is implemented 
employing the input parameters extracted from the EH 
calculations. The results are shown in Table [T] and Figs. 
M and El 

First, let us focus on the prototypical three-block sym- 
metrical structure PDA52/PA65/PDA52. The potential 
profiles are similar to the ones shown in Fig. [3] The 
well (PA) region and the entire molecule have lengths 
L w = 16.02 nm and L m = 67.05 nm, respectively. The 
number of bound states, N, for a particle of mass m con- 
fined in a symmetrical finite well of depth V satisfies the 

inequality N — 1 < (2mV ^ /2L '" < iVP Employing this 
formula, with the average of the hole (electron) effective 
masses of PA and PDA, the EMO method estimates the 
appearance of four (three) states confined within the VB 
(CB) potential well. This is precisely what the numer- 
ical solution of Eq. (JT9J) yielded. In Table [i] the ener- 
gies of these states are compared with the EH MO en- 
ergies. The agreement of the EMO hole (electron) en- 
ergies with the ones of the corresponding EH HOMO-3, 
HOMO-2, HOMO-1 and HOMO (LUMO, LUMO+1 and 
LUMO+2) is very good, although it increasingly deteri- 
orates as this variable decreases (increases). Table [I] also 
compares the EMO energy of an electron state located 
slightly above the intramolecular barrier with the EH 
LUMO+3 energy. A very good agreement is found here 
too, which, interestingly, is better than for LUMO+1 and 
LUMO+2. The differences are of the order of 1 — 10 meV, 
which is very encouraging. 

Since the EMO's and EH MO's are one-dimensional 
and three-dimensional functions, respectively, they can- 
not be directly compared. To devise an useful comparison 
criterion, the shape of an MO along the longitudinal di- 
rection of the molecule is described in terms of the set of 
parameters 



4" 



(22) 



where c$ is the coefficient of the r-th atomic basis or- 
bital centered at the j-th carbon core in the q-th EH MO. 
The contributions of the Is atomic orbitals centered at 
the hydrogen nuclei are not taken into account. The set 
of atomic indices j and the parameter can be in- 
terpreted as a discretization of the z coordinate and the 
probability density associated with the q-th MO at the 
j-th position, respectively. 

Fig. [5] shows a comparison of such MO longitudinal 
shapes (appropriately scaled) with the EMO probabil- 



ity densities. It is seen, firstly, that the EH MO's display 
long-wavelength large-amplitude envelopes, accompanied 
by short-wavelength low-amplitude oscillations that can 
be attributed to the band-edge background Bloch func- 
tions, as Eqs. (14) and (15) predict. Secondly, the EH 



MO's confined within the well (PA) region exhibit the 
same behavior as the square- well eigenfunctions. Thirdly, 
the EMO's describe very precisely the long-range varia- 
tions of the EH MO's, except for deviations of HOMO-3, 
LUMO+1 and LUMO+2 within the barrier regions, and 
LUMO+3 within the well region. 

Let us assess qualitatively the contribution of each 
one of the main approximations involved in the theo- 
retical development to the origins of the observed differ- 
ences in energy and spatial behavior between the EMO's 
and EH MO's. In the first place, the fact that no ap- 
preciable spatial deviations are found in the immedi- 
ate neighborhoods of the hetero junctions indicates that 
the uncoupled approximation and the artificial bound- 
ary conditions work very well for this system. This is 
likely due to the structural similarity of PA and PDA. 
In the second place, the observed agreement between the 
EMO and EH MO wavefunctions within the classically- 
allowed regions and the deviations thereof within the 
classically-forbidden regions corroborate the line of argu- 
mentation preceding Eqs. (14) and (15), concerning the 
single-envelope approximation. In particular, an EMO 
is more delocalized within the classically-forbidden re- 
gions than the corresponding EH MO because there the 
latter is more accurately represented by the convolution 
(z) with the Fourier transform to 



^(Rpda) 



product of F, 

z-space of u^ DA \k,r) (see Eqs. Q and |ll])), than 

by the algebraic product F^ PDA \z)un DA \k, f) (Eqs. 
( 14 ) and ( 15 )). Such convolution yields a narrower wave- 
function because Un DA \k, r) is confined to the first BZ 
in fc-space, so its Fourier transform to z-space is local- 
ized within the length of the PDA unit cell, i a . Nat- 
urally, these deviations are more noticeable for states 
close to the edges of the wells than for deeply-lying ones, 
since the wavefunctions of the former penetrate more 
into the classically-forbidden regions. In the third place, 
the hole and electron MO's confined inside the well have 
mostly PA character and their energies are seen to lie 
inside the parabolic regions of the VB and CB disper- 
sions of PA, respectively (see Table [I] and Figs. [2] [3] and 
[4j. Hence, the error introduced by the parabolic EMA 
is very small for these states. On the other hand, the 
LUMO+3 has a dominant PDA character because it is 
not confined within the well. Since its energy lies slightly 
above the barrier, it lies deeply within the parabolic re- 
gion of the CB of PDA. The very slight deviations ob- 
served across the well region must be due to small errors 
in the parabolic EMA for PA, since this state is relatively 
far from the bottom of the well. This explains why the 
agreement for this state is better than for LUMO+1 and 
LUMO+2. 

Second, let us consider the shorter three-block sym- 
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FIG. 5. Comparison of the longitudinal MO shapes (gray wiggly lines) and the EMO probability densities (red smooth lines) 
for (a) LUMO and HOMO, (b) LUMO+1 and HOMO-1, (c) LUMO+2 and HOMO-2, and (d) LUMO+3 and HOMO-3 of 
PDA52 /PA65 /PDA52 (see Table |T]| . Functions for holes are displayed "upside down" . The vertical dashed lines indicate the 
positions of the intramolecular barriers. The horizontal scale is in rim. 



metrical structure PDA 26 /PA4o/PDA 26 . The well (PA) 
region and the entire molecule have lengths L w — 9.86 
nm and L m = 35.37 nm, respectively. Employing the 
same considerations as for the previous QD, the EMO 
method estimates the appearance of two states confined 
in each potential well. Again, this is precisely what the 
numerical solution of Eq. ( 19 ) produced. Table [I] presents 
the energies of these states plus the energies of one addi- 
tional hole state and one additional electron state located 
below and above the corresponding barriers. Compari- 
son with the energies of the corresponding EH HOMO- 
2, HOMO-1, HOMO, LUMO, LUMO+1 and LUMO+2, 
shows qualitatively the same behavior as for the previ- 
ous QD, although now the deviations tend to be about 
twice as large. This is expected, since now the potential 
well and the molecule are shorter, causing the hole and 
electron states to be pushed closer to the edges of their re- 
spective wells, and making the influence of the interfaces 
more noticeable. Fig [6] displays the comparison between 
the EMO probability densities and the EH MO shapes. 
Again, the qualitative behaviors observed are analogous 
to the ones of the previous QD. Quantitatively, it can 
be seen that the deviation of the (unconfined) LUMO+2 
across the well (PA) region is larger than the one of the 
LUMO+3 of the previous case. This is because the latter 
is closer to the LUCO of PA, so the parabolic EMA in this 
region works better for it. It is also interesting to notice 
that the EH HOMO-2 is bound with respect to the well, 



while the corresponding EMO is barely unbound, lying 
only 2 meV above the barrier. This effect is attributed 
to a small deviation from the parabolic approximation in 
the effective mass. 

Third, let us consider the two-block asymmetrical 
structure PDA 2 o/PA 2 o, which, according to Fig. Q, 
exhibits a well within the PA region. The well and 
the entire molecule have lengths L w — 4.92 mn and 
L m = 14.72 nm, respectively. Table [l] shows that now 
only the HOMO and LUMO are confined within the well 
(PA) region. The degree of agreement between the EH 
and EMO energies for this structure is slightly lower than 
for PDA 2 6/PA4o/PDA 2 6, due to its shorter length, but 
is still very good. The EMO probability densities and 
the EH MO shapes (not shown) display the expected de- 
formations associated with the asymmetry of the struc- 
ture, and qualitatively agree to the same degree as in 
PDA 26 /PA 40 /PDA 26 . 

Finally, let us consider the one-block (single 
oligomer) structure PAgs, which can be visual- 
ized as the "heterostructures" vacuum/PAgs /vacuum, 
vacuum/PA 2 o/PA 2 5/PA 2 o /vacuum or any other ap- 
propriate partition of 65 PA monomers. A com- 
parison of this system with the true heterostructure 
PDA5 2 /PA65 /PDA5 2 helps to further assess the effects of 
the uncoupled approximation and the artificial boundary 
conditions at oligomer-oligomer and oligomer- vacuum in- 
terfaces. In Table U it is first observed that the electron 



Envelope molecular-orbital theory of extended systems I 



12 




FIG. 6. Analogous to Fig. [5] for (a) LUMO and HOMO, (b) LUMO+1 and HOMO-1, and (c) LUMO+2 and HOMO-2 of 
PDA 26 /PA 4 o/PDA26. 



and hole levels of PA 6 s have been slightly pushed up 
and down, respectively, in comparison with the corre- 
sponding confined levels of PDA52/PA65/PDA52, which 
is expected because the PA-vacuum barriers are much 
larger than the PA-PDA barriers, causing a stronger con- 
finement. Second, the degree of agreement between the 
EMO and the EH MO energies is somewhat lower than 
for PDA 52 /PA 65 /PDA52, which makes sense because the 
potential and effective-mass jumps are much larger at the 
PA-vacuum interface than at the PA-PDA interface, so 
the effects of the uncoupled approximation and the arti- 
ficial boundary conditions should be more noticeable. 



B. DFT Input Parameters 

In this subsection the EMO method is implemented for 
the smaller structures, employing the input parameters 
extracted from the DFT calculations. The results are 
shown in Table [TI) In this case, a comparison between the 
EMO probability densities and the Kohn-Sham orbital 
shapes will not be provided. 

First, let us consider again the two-block asymmetri- 
cal structure PDA 2 o/PA 2 o, which, according to the DFT 
alignments (Fig. now exhibits a well for electrons 
within the PDA region and a well for holes within the 
PA region. Table |TT] shows that according to the EMO 
method there are now two confined hole states within the 
well (PA) region, while with the EH parameters only one 
had been found. This comes about because the DFT off- 
set is larger than the EH one. For electrons, three states 
confined within the well (PDA) region are found. The 
degree of agreement between the DFT and EMO ener- 
gies is about the same for some states and lower for other 
states in comparison with the one between the EH and 



EMO energies of the previous subsection. 

Second, let us consider the one-block (single 
oligomer) structure PDA20, which can be visual- 
ized as the "heterostructures" vacuum/PDA2o /vacuum, 
vacuum/PDA5/PDA 10 /PDA 5 /vacuum or any other ap- 
propriate partition of 20 PDA monomers. It is seen that 
the level of agreement between the DFT and EMO en- 
ergies is slightly lower than for PDA2o/PA2o- This is 
expected for reasons analogous to the ones argued for 
the last structure studied in the previous subsection. 
Furthermore, it is expected that the PDA20 electron 
states should be slightly pushed up in comparison with 
the electron states confined within the PDA region in 
PDA20/PA20, because the PDA- vacuum barrier is much 
larger than the PDA-PA one, causing a stronger confine- 
ment. It is observed that the EMO results display the 
correct trend of behavior, whereas the DFT ones do not. 
This comes as no surprise, since DFT/PBE is known 
to be rather inaccurate for conjugated systems.^ (An 
analogous comparison for holes is not possible because in 
PDA20/PA20, these are confined within the PA region.) 

This section is closed with the following additional 
comments: First, the results of calculations on smaller 
structures (not shown) exhibited the same trends of be- 
havior, in particular the expected deterioration with de- 
creasing oligomer size of the agreement between the EMO 
and atomistic results. On the other hand, it should be 
evident by now that results for very large, mesoscopic 
heterostructures are likely to be equally, or even more, 
accurate, than the ones shown here. Second, no local- 
ized states at oligomer-oligomer or oligomer-vacuum in- 
terfaces were found in these structures from the atomistic 
calculations within the energy region examined, although 
it is pos sible that such states appear at higher or lower 
energies! 34 ^ 35 ! The uncoupled approximation and the ar- 
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tificial boundary conditions used in this work preclude a 
prediction of this kind of state. 



VI. CONCLUSIONS AND PERSPECTIVES 

The EMO theory introduced in this paper, which 
constitutes an extension for (finite) molecules of the 
envelope-function approximation widely employed in 
solid-state physics, is seen to provide an attractive con- 
ceptual and computational framework for the description 
of the course-grained (non-atomistic) electronic structure 
of extended quasi-linear heterostructures. Its simplicity, 
as compared with atomistic electronic-structure meth- 
ods, stems from the mapping of the three-dimensional 
single-particle eigenvalue equations into one-dimensional 
Schrodinger-like equations, whose eigenf unctions de- 
scribe automatically the envelopes of the MO's. The 
practical implementation of this approach requires four 
input parameters for each parent polymer, namely the 
HOCO (top of valence band) and LUCO (bottom of con- 
duction band) energies and the hole and electron effective 
masses (in some cases, the LOCO energies of the lateral 
parent polymers may also be needed), which, ideally, can 
be extracted from experimental measurements or, in lieu 
of these, from electronic bandstructure calculations. It 
should be born in mind that the need for these parame- 
ters arises not because the EMO approach is intrinsically 
a semiempirical model, but because it is an effective the- 
ory. 

As an illustration, the approach was applied to model 
7r-conjugated co-oligomers constituting so-called quan- 
tum dots, employing input parameters extracted from 
EH and DFT bandstructure calculations. To meaning- 
fully assess the internal consistency and accuracy of the 
EMO approach, its predictions were compared with the 
results of concomitant atomistic molecular electronic- 
structure calculations. For high-lying holes and low-lying 
electrons, the agreement between the EMO and atom- 
istic energies was very good (within the order of 1 — 10 
meV). In addition, the envelope MO's were found to 
describe very precisely the long-wavelength oscillations 
of the full MO's, except when these penetrate deeply 
into classically-forbidden zones, where small deviations 
are observed. The origins of all these discrepancies were 
qualitatively accounted for. 

The approximations involved in the formal develop- 
ment were introduced on heuristic grounds and validated 
a posteriori. The main shortcoming resulting from these 
approximations is the incapability of predicting states lo- 
calized at the heterojunctions. A quantitative estimation 
of errors and the prediction of such localized states re- 
quire the examination of the couplings between bands 
and the reconstruction effects at the interfaces. These 
tasks will be undertaken in future reports. 

It is of considerable interest to ext end th e EMO the- 
ory for the calculation of the optical 4 ^ 6 ^* anc | charge- 
transport!^ properties of conjugated molecular struc- 



tures. Work is under way in these directions! 61 * 62 ! 

It would be useful also to generalize the EMO method 
to systems with diverse compositions and higher dimen- 
sionalities, for example graphene-based str uctur e d 6 * 63 * 64 ! 
and hyper-branched, or dendritic, molecules.* 4 ^ 

As far as applications are concerned, the ability to de- 
sign molecular structures with properties specifically en- 
gineered for the construction of bulk or molecular elec- 
tronic and optoelectronic devices is very important! 46 * 66 * 
The EMO method may constitute a useful aid for 
this task, since it readily permits the tailoring of 
the energy levels and wavefunctions, by controlling 
the sizes of the constituent fragments, avoiding the 
costly and tedious performance of extensive electronic- 
structure calculations every time. For example, it 
has already been applied to the design of nanomet- 
ric molecular architectures^ for quantum computing^!, 
molecular sup erlattic es^ with potential applications in 
optoelectronics^ 4 ^ 6 -* and nanotube-based double-barrier 
heterostructures for chemical sensing.^. Moreover, it is 
desirable to develop methodologies for the design of par- 
ent materials with bandgaps and effective masses tailored 
for specific applications; progress in this direction is also 
under way! 6 ^* 

There is an obvious similarity between the EMO the- 
ory, if applied to a single oligomer, and the FEMO 
mo del E2H2S1 Clearly, in the former the effective box-like 
potentials and masses for electrons and holes arise as 
a matter of course from a first-principles formalism, 
whereas in the latter the box potential and effective mass 
for electrons are introduced on grounds of intuitive, al- 
beit very clever, considerations. This lack of a rigor- 
ous derivation of the FEMO model hampers its system- 
atic improvement. Nevertheless, now it appears that this 
model can, in fact, be easily deduced from the EMO the- 
ory by making a few additional approximations. 
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